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Abstract Distortion risk measures are extensively used in finance and insurance ap¬ 
plications because of their appealing properties. We present three methods to construct 
new class of distortion functions and measures. The approach involves the composting 
methods, the mixing methods and the approach that based on the theory of copula. 
Subadditivity is an important property when aggregating risks in order to preserve the 
benefits of diversification. However, Value at risk (VaR), as the most well-known exam¬ 
ple of distortion risk measure is not always globally subadditive, except of elliptically 
distributed risks. In this paper, instead of study subadditivity we investigate the tail 
subadditivity for VaR and other distortion risk measures. In particular, we demonstrate 
that VaR is tail subadditive for the case where the support of risk is bounded. Various 
examples are also presented to illustrate the results. 
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1 Introduction 


A risk measure p is a mapping from the set of random variables A, standing for risky 
portfolios of assets and/or liabilities, to the real line M. In the subsequent discussion, 
positive values of elements of X will be considered to represent losses, while negative 
values will represent gains. Distortion risk measures are a particular and most important 
family of risk measures that have been extensively used in finance and insurance as capital 
requirement and principles of premium calculation for the regulator and supervisor. Sev¬ 
eral popular risk measures belong to the family of distortion risk measures. For example, 
the value-at-risk (VaR), the tail value-at-risk (TVaR) and the Wang distortion measure. 
Distortion risk measures satisfy a set of properties including positive homogeneity, trans¬ 
lation invariance and monotonicity. When the associated distortion function is concave, 
the distortion risk measure is also subadditive (Denneberg, 1994; Wang and Dhaene, 
1998; Wirch and Hardy, 2001). VaR is one of the most popular risk measures used in risk 
management and banking supervision due to its computational simplicity and for some 
regularity reasons, despite has some shortcomings as a risk measure. For example, VaR 
is not a subadditive risk measure (see, for instance, Artzner et al. (1999), Denuit et al., 
(2006)), it only concerns about the frequency of risk, but not the size of risk. TVaR, 
although being coherent, concerns only losses exceeding the VaR and ignores useful infor¬ 
mation of the loss distribution below VaR. Clearly, it is difficult to believe that a unique 
risk measure could capture all characteristics of risk, so that an ideal measure does not 
exist. Moreover, since risk measures associate a single number to a risk, as a matter of 
fact, they cannot exhaustively all the information of a risk. However, it is reasonable to 
search for risk measures which are ideal for the particular problem under investigation. 
As all the proposed risk measures have drawbacks and limited applications, the selection 
of the appropriate risk measures continues to be a hot topic in risk management. 

Zhu and hi (2012) introduced and studied the tail distortion risk measure which was 
reformulated by Yang (2012) as follows. For a distortion function g, the tail distortion 
risk measure at level p of a loss variable X is defined as the distortion risk measure with 
distortion function 



Some properties and applications can be found in Mao, Lv and Hu (2012), Mao and Hu 
(2012) and Lv, Pan and Hu (2013). 
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As an extension of VaR and TVaR, Belles-Sampera et al. (2014a) proposed a class 
of new distortion risk measures called Glue VaR risk measures, which can be expressed as 
a combination of VaR and TVaR measures at different probability levels. They obtain 
the analytical closed-form expressions for the most frequently used distribution functions 
in financial and insurance applications, while a subfamily of these risk measures has 
been shown to satisfy the tail-subadditivity property which means that the benefits of 
diversification can be preserved, at least they hold in extreme cases. The applications 
of Glue VaR risk measures in capital allocation can be found in the recent paper Belles- 
Sampera et al. (2014b). 

Cherubini and Mulinacci (2014) propose a class of distortion measures based on con¬ 
tagion from an external “scenario” variable. The dependence between the scenario and 
the variable whose risk is measured is modeled with a copula function with horizontal 
concave sections, they give conditions to ensure that coherence requirements be met, and 
propose examples of measures in this class based on copula functions. 

The first purpose of this paper is to construct new risk measures following Zhu and Li 
(2012), Belles-Sampera et al. (2014a) and Cherubini and Mulinacci (2014). The newly 
introduced risk measures are included the tail distortion risk measure and the GlueVaR as 
specials. The second goal of the paper is to investigate the tail asymptotics of distortion 
risk measures for the sum of possibly dependent risks with emphasis on VaR. The rest of 
the paper is organized as follows. We review some basic definitions and notations such as 
distorted functions, distorted expectations and distortion risk measures in Section 2. In 
Section 3 several new distortion functions and risk measures are introduced. In Section 
4 we investigate the tail asymptotics as well as subadditivity/superadditivity of VaR. 
Finally, in Section 5 we analyze the subadditivity properties of a class of distortion risk 
measures. 


2 Distortion risk measures 

2.1 Distorted functions 

A distortion function is a non-decreasing function g : [0,1] —> [0,1] such that p(0) = 
0,g(l) = 1. Since Yaari (1987) introduced distortion function in dual theory of choice 
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under risk, many different distortions g have been proposed in the literature. Here we 
list some commonly used distortion functions. A summary of other proposed distortion 
functions can be found in Denuit et al. (2006). 

• g(x) — l( x> i_ p )(a;), where the notation 1 a to denote the indicator function, which 
equals 1 when A holds true and 0 otherwise. 

• g(x) = min{y^, 1}. 

• Incomplete beta function g(x) = ] 3 (tb) fo ta 1 (1 — t) b 1 dt 1 where a > 0 and b > 0 

are parameters and /3(a, b ) = / 0 ' t a_1 (l — Setting 6 = 1 gives the power distortion 

g{x) = x a . Setting a = 1 gives the dual-power distortion g(x) = 1 — (1 — x) b . 

• The Wang distortion g(x) = < h( < h _1 (x) + $ _1 (p)),0 < p < 1, where $ is the distri¬ 
bution function of the standard normal. 

• The lookback distortion g(x) = x p ( 1 — p In x),p G (0,1]. 

Obviously, every concave distortion function is continuous on the interval (0,1] and 
can have jumps in 0. In contrast, every convex distortion function is continuous on the 
interval [0,1) and can have jumps in 1. For a distortion function g, if there exists a to > 0 
such that g(t 0 ) = 0, then g is not concave; if there exists a t\ < 1 such that g(t\) = 1, 
then g is not convex. The identity function is the smallest concave distortion function 
and also the largest convex distortion function; go(x) := l(a;>o) is concave on [0,1] and 
is the largest distortion function. g°(x) := l( x= i) is convex on [0,1] and is the smallest 
distortion function. For 0 < p < 1, we remark that g\(x) := min-fy^, 1} is the smallest 
concave distortion function such that gi(x) > l( 3 ,>i_ p )(x). In fact, we consider a concave 
distortion function g such that g(x) > l( a;> i_p), then g = 1 on (1 —p, 1], As g is concave, 
it follows that g(x) > for x < 1 — p, and thus g{x) > min {, 1} for 0 < x < 1. 
Any concave distortion function g gives more weight to the tail than the identity function 
g(x) = x, whereas any convex distortion function g gives less weight to the tail than the 
identity function g{x) = x. 
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2.2 Distorted risk measures 


Let (0, F, P) be a probability space on which all random variables involved are defined. 
Let Fx be the cumulative distribution function of random variable X and the decumula- 
tive distribution function is denoted by F x , i.e. F x (x ) = 1 — F x (x) = P(X > x). Let g 
be a distortion function. The distorted expectation of the random variable X, notation 
p g [X], is defined as 

/*+oo pO 

Pg[X] = / g(F x (x))dx+ / \g(F x (x)) - l]dx, 

J 0 J — oo 


provided at least one of the two integrals above is finite. If X a non-negative random 
variable, then p g reduces to 


r»+oo 


Pg[X] = / g(F x (x))dx. 


Jo 

From a mathematical point of view, a distortion expectation is the Choquet integral 
(see Denneberg (1994)) with respect to the nonadditive measure p = g o P. That is 
p g [X\ = f Xdp. In view of Dhaene et al. (2012, Theorems 4 and 6) we know that, when 
the distortion function g is right continuous on [0,1), then p g [X] may be rewritten as 


p,[X] = / Va Rt_ t [X)ig(q), 

•'[Ml 


where VaR + p[X] = sup{x|Fx(x) < p}, and when the distortion function g is left contin¬ 
uous on (0,1], then p g [X] may be rewritten as 

p g [X\ = [ V aRi- q [X]dg (q) = I VaR q [X]dg(q), 

J[ 0,1] J[ 0,1] 

where VaR p [X] = inf{x|Fx(x) > p} and g(q) := 1 — g( 1 — q) is the dual distortion of g. 
Obviously, g — g, g is left continuous if and only if g is right continuous; g is concave if 
and only if g is convex. The distorted expectation p g [X] is called a distortion risk measure 
with distortion function g. Distortion risk measures are a particular class of risk measures 
which as premium principles were introduced by Deneberg (1994) and further developed 
by Wang (1996, 2000) among others. As it is well known, the mathematical expectation, 
E[X], is a distortion risk measure whose distortion function is the identity function. If g 
is concave, then 


r»+00 


Pg[X\ > / F x {x)dx + / [F x (x) - 1 }dx = E[X], 
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and if g is convex, then 


/»+oo /*0 

p g [X\ < / F x (x)dx + / [F x (x) - 1 \dx = E[X}. 

J 0 J — oo 

Distortion risk measures satisfy a set of properties including positive homogeneity, 
translation invariance and monotonicity. Hardy and Wirch (2001) have shown that a risk 
measure based on a distortion function is coherent if and only if the distortion function 
is concave. A risk measure is said to be coherent if it satisfies the following set of four 
properties (see, e.g., Arztner et al. 1997 and 1999): 

(M) Monotonicity: p(X) < p(Y ) provided that P(X < Y) = 1. 

(P) Positive homogeneity: For any positive constant c > 0 and loss X,p(cX) = cp(X). 

(S) Subadditivity: For any losses A", Y, then p(X + Y) < p(X) + p(Y). 

(T) Translation invariance: If c is a constant, then p(X + c) = p(X) + c. 

It is furthermore shown by Artzner et al. (1999) that all mappings satisfying the above 
properties allow a representation: 


p(X) = supBJX], 

peV 

where V is a collection of ‘generalised scenarios’. A risk measure p is called a convex risk 
measure if it satisfies monotonicity, translation invariance and the following convexity 
(C): 

p(XX + (1 - A)Y) < Xp(X) + (1 - A)p(Y), 0 < A < 1. 

Clearly, under the assumption of positive homogeneity, monotonicity and translation in¬ 
variance, the convexity of a risk measure is equivalent subadditivity. 


The most well-known examples of distortion risk measures are the above-mentioned 
VaR and TVaR, corresponding to the distortion functions, respectively, are g{x) = l( a:> i-p) 
and g{x) = min j, l|. Notice that TVaR p [A] can be alternatively expressed as the 
weighted average of VaR and losses exceeding VaR: 

TVaR p [A] = VaR p [X] + 1 ~ F x( Va R P [X]) E ^ _ VaR ^ x ^ x > Va R p [X}\. (2.1) 


For continuous distributions, TVaR coincide with the expected loss exceeding p-Value-at 
Risk, i.e., the mean of the worst (1— p)100% losses in a specified time period which defined 
by 


CTE P [X] = E[X\X> VaR p [X]\. 
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If X is a real valued random variable and 0 < p < 1, then we say that q is an p-quantile 
if P[X < q] < p < P[X < q}. By definition, VaR p [X] is the lower p-quantile of the 
r.v. X and VaR p [X] is the upper p-quantile of the r.v. X. VaR p [X] is a left-continuous 
nondecreasing function having VaRo[X] as the essential inhmum of X , possibly — oo, 
VaR p [X] is a right-continuous nondecreasing function having VaRi[X] as the essential 
supermum of X, possibly +oo. It is easy to see that VaR p [X] < VaR p [X], there are 
at most countably many values of p G [0,1] where VaR p [X] and VaR p [X] differ (see, 
Dhaene et al. (2012)). Moreover, VaR p [X] = VaR p [X] if, and only if Fx(x) = p for 
at most one x, which equivalent to Fx(') is strictly increasing. The risk measures VaR 
and VaR + satisfy axioms (M), (P), and (T), but not (S) and (C) (except in some special 
cases, such as in the multivariate normal distributions or more generally multivariate 
elliptical distributions), and hence is not coherent in the sense of Artzner et al. (1999). 
Despite suffers from some serious limitations, VaR is still the standard of industry and 
regulatory for the calculation of risk capital in banking and insurance. For example, the 
Basel Committee on Banking Supervision introduced a 99% Value at Risk requirement, 
based on a 10-day trading horizon. The TVaR improves the VaR as a measure of risk by 
also taking into account the magnitude of loss beyond the VaR. That is TVaR measures 
average losses in the most adverse cases rather than just the minimum loss, as the VaR 
does. Therefore, risk assessment based on the TVaR have to be considerably higher than 
those based on VaR. The importance of TVaR is also seen from a result of Kusuoka 
(2001), who proved that TVaR p is the smallest law invariant coherent risk measure that 
dominates VaR p . Unlike VaR, the distortion function associated to the TVaR is concave 
and, then, the TVaR is a coherent risk measure in the sense of Artzner et al. (1999). It 
means that TVaR is a subadditive risk measure (see, for instance, Denuit et ah, 2006). 
In the literature, the TVaR is sometimes called the expected shortfall. Although TVaR 
is one of the best coherent risk measures, however, TVaR reflects only the mean size of 
losses exceeding the VaR. It ignores the useful information in a large part of the loss 
distribution, and consequently lacks incentive for mitigating losses below the quartile 
VaR. Moreover, it does not properly adjust for extreme low-frequency and high-severity 
losses, since it only accounts for the mean value (not higher moments). A recent paper by 
Frittelli et al. (2014) has proposed a new risk measure, the lambda value at risk AVaR) as 
a generalization of the VaR. The novelty of the AVaR lies in the fact that the confidence 
level can change and adjust according to the risk factor profit and loss. 

Detailed studies of distortion risk measures and their relation with orderings of risk 
and the concept of comonotonicity can be found in, for example, Wang (1996), Wang and 
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Young (1998), Hiirlimann (1998), Hua and Joe (2012) and the references therein. The 
following lemma will be used in proofs of later results, which characterizes an ordering of 
distortion risk measures in terms of their distortion functions. 

Lemma 2.1. (Belles-Sampera et al. (2014b)). If g(x) < g*(x ) for x G [0,1], then 
p g [X] < p g *[X] for any random variable X. 


3 Generating new distortion functions and measures 


Distortion functions can be considered as a starting point for constructing families of dis¬ 
tortion risk measures. Thus, constructions of distortion functions play an important role 
in producing various families of risk measures. Using the technique of mixing, composition 
and copula allow the construction of new class of distortion functions and measures. 


3.1 Composting methods 

The first approach to construct distortion functions is the composition of distortion func¬ 
tions. 

Let be distortion functions, define f\(x) = h\{x) and composite functions 

fn(x) = fn-i(h n (x)), n = 1,2,---. ft is easy to check that f n (x),n = 1,2, ••• are all 
distortion functions. If h 1 ,h 2 ,--- are concave distortion functions, then each f n (x) is 
concave and satisfies that 

,/ l < f'l < ,/3 < • • • 

and 

lim f n (x) = 1( X>0 ), x G [0,1]. 

71—>• OO 

The associated risk measure satisfies (by Lemma 2.1) 

Ph\X)<Pt,[X}<P h [X}<--- 

and 

lim Pf n [X] = VaRi[X] = esssup(A"). 

71—>• OO 

If hi,h 2 ,--- are convex distortion functions, then each f n (x) is convex and satisfies that 

fi > h > h > • • • 
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and 


lim f n (x) = l(*=i), x e [0,1]. 

n—>-oo 

The associated risk measure satisfies (by Lemma 2.1) 

piAx\> Ph\x)> piA *]>••• 

and 

lim Pf n [X) = VaRo[X] = essinf(Y). 

n—> oo 

Consider two distortion functions gi and g 2 . If 

, f if 0 < x < 1-p, 

| 1, if 1 — p < x < 1, 

then we get 

f gi ( t 2 -') , if 0 < x < 1 — p, 
g P (x) ■■= gMx)) = { y v 1 ^’ - - 

1, if 1 — p < x < 1. 

The corresponding risk measure p ;]p [X is the tail distortion risk measure which was first 
introduced by Zhu and Li (2012), and was reformulated by Yang (2012). In particular, 
on the space of continuous loss random variables X, 

POO 

Pg p [X] = g p (1 - P(X < x\X > VaR p [X})) dx. 


If g\(x) = x r , 0 < r < 1 and 


/ n i i—p’ if 0<x< 1-p, 

9 ^ x ) = { 1 P ^ . . 

1, it 1 — p < x < 1, 


then 


gi2 (x) := gi(g 2 (x)) = 


, if 0 <x< 1-p, 
, if 1 — p < x < 1, 


and 


( \ J T=C lf 0 - X “P) r > 

g 2 iyx) :=g 2 {g 1 {x)) = { \ p 

1, if (1 -p)~ < X < 1. 

Clearly, g\ < g 2 \ and g 2 < g \ 2 , so that, by Lemma 2.1, p gi [X] < p 921 [X] and p 92 [X] < 
Pgi2 [-^1 • 
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In practice, sometimes one needs distort the initial distribution more than one times. 


Example 3.1 Consider two risks X and Y with distributions, respectively, are: 


F x {x) = < 


and 


F y (x) = 


0, 

if x < 0, 

0.6, 

if 0 < x < 100, 

0.975, 

if 100 < x < 500, 

1, 

if x > 500, 

o, 

if x < 0, 

0.6, 

if 0 < x < 100, 

0.99, 

if 100 < x < 1100, 

1, 

if x > 1100. 

—VaR 0 .96[X] —100, VaRo.gsfY] 


TVaR can be calculated by formula (2.1): 


TVaRo.95[X]=TVaRo.95[Y]=300, TVaRo.ge[X]=TVaRo .96 [Y]=350. So that when a = 0.95 
and /? = 0.96, according to the measures of VaR and TVaR, both X and Y bear the same 
risk! However, the maximal loss for Y (1100) is more than double than for loss X (500), 
clearly, risk Y is more risky than risk X. Now we consider distortion expectation p 9p with 


9i(x) = 92 (x) 


if 0 < x < 1 — p, 
1, if 1 — p < x < 1. 


One can easily find that, with p = 0.95, p 9p [X] = 500 and p 9p [Y] = 1100. 


3.2 Mixing methods 

One of the easiest ways to generate distortion functions is to use the method of mixing 
along with finitely distortion functions or infinitely many distortion functions. Specifically, 
if g w (w G< a,b >) is a one-parameter family of distortion functions, ^ is an increasing 
function on < a,b > such that f b> d'i/j(w) = 1, then the function g = f <a h> g w dtp(w) is 
a distortion function, the associated risk measure is given by 

pAX\ = /' x\dpw). (3.1) 

J <a,b> 
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In particular, if if is discrete distribution, then (3.1) can be written as the form of convex 
linear combination g = w^i (Wi > 0, JA w i — 1) , the associated risk measure is given 
by 

Pg[ x ] = Y^ w iPgi[ x \- ( 3 - 2 ) 

i 


The following lemma is well known (cf. Kriele and Wolf (2014, Theorem 2.1, P.33)). 

Lemma 3.1. If all p 9w (w E< a,b>) are monotone, positively homogeneous, subadditive 
and translation invariant, then p g [X] also has the corresponding properties. That is, if all 
g w (w G< a,b>) are coherent, then p g [X] is also coherent. 


Now we list three interesting special cases: 


• If [a, b) = [0, oo), gi(x) = 1 — (1 — x)\i > 1 and Wi > 0 = then p g [X\ in 

(3.2) is coherent since gi(x) is concave. As in Tsukahara (2009), if we take w, from Bin 
(1, 6) (0 < 9 < 1), then g g {u ) = u + uO — u 2 6. If we take 


Wi = 


y, 0 > 0, 


then 


(«) = 


l)i\ 


e e (l — e 9u ) 


Also, if take Wi = (1 — 0) 1 l 8 (0 < 0 < 1), the geometric distribution, then 


g e (u) = 


u 


u + 6(1 — u ) ’ 

which is the proportional odds distortion; see Example 2.1 in Cherubini and Mulinacci 
(2014). 


If [a, b] = [0,1], p 9w = VaR w [X] and dtp(w) = 4>(w)dw, then p g [X\ in (3.1) reduces 


to 


P 4 >[X] = / VaR w [X\cj)(w)dw, 


(3.3) 


which is spectral risk measure (see Acerbi 2002, 2004). Here <f is called a weighting 
function satisfies the following properties: 0 > 0, (f(w)dw = 1. The following lemma 
gives a sufficient condition for p^X] to be a coherent risk measure (cf. Kriele and Wolf 
(2014)). 


Lemma 3.2. Spectral risk measure p ( j > [X\ is coherent if cf is (almost everywhere) mono¬ 
tone increasing. 


11 





Clearly, there exists a one-to-one correspondence between distortion function g and 
weighting function q i>, namely, g(l — t) = 1 — f* (f)(s)ds. Obviously, g is concave if, and 
only if 0 is (almost everywhere) monotone increasing. Two well-known members of this 
class are the VaR and the TVaR. The associated weight functions are 0(iu) = S p (w) and 
-p-X (w > p ), respectively. Here S p (w) is a Dirac delta function that gives the outcome 
a — p an infinite weight and gives every other outcome a weight of zero. From Lemma 
3.2, TVaR is coherent since = -^\(w > p) is monotone increasing. By contrast, 
(j>(w) = S p (w) is not monotone increasing, hence VaR is not coherent. Both of these 
measures use only the tail of the distribution. 


• If [a, b] = [0,1], p 9w = TVaR w [X] and 0 = /i is a probability measure on [0,1], then 
p g [X\ in (3.1) reduces to 

p p [X\= I TVaR w [X}dp(w ), (3.4) 

Jo 

which is the weighted TVaR (see Cherny (2006)). TVaR p is a special weighted TVaR 
with p(w) = l(tc > p). According to Lemma 3.1, since each TVaR^AT] is coherent 
risk measure, the weighted TVaR is coherent risk measure. The weighted TVaR can be 
rewritten as the form of spectral risk measure as following: 


p p [X] = / TVaR w [X}dp(w) = 


VaR q [X] 


1 — w 


VaR q [X]dq ) dp(w 


1 — w 


dp(w) ) dq (by the Fubini theorem) 


= / VaR q [X](j>(q)dq 

Jo 

= [ VaRi_ q [X]dg(q), 

Jo 

where, g is a function with g( 0) = 0 and satisfies 

g\ 1 -q) = <K<?) = 


9 1 


1 — w 


dp{w). 


Because 0(g) is increasing function of g, it follows from Lemma 3.2 that the weighted 
TVaR p p [X] is coherent. Or, equivalently, g\q) is decreasing function of g, i.e. g is a 
concave function, moreover, g is increasing and 


0 (!) = / g'( 1 -w)dw = dq 

Jo Jo 

r l i r l 

dp(w) / dq 


1 — w 


dp(w ) 


1 — w 


dp(w) = 1. 
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so that g is a concave distortion function, and hence the weighted TVaR p^[X\ is coherent. 


Conversely, the distortion measure with concave distortion function g can be expressed 
by the weighted TVaR. In fact, note that (j)(q ) = g'( 1 — q ) is monotone increasing, we 
define a measure z/([0 ,q]) = 4>(q). As in the proof of Theorem 2.4 in Kriele and Wolf 
(2014) we have 


p 9 [X\ = - VaR w [X]dg(l-w) 

Jo 

= [ VaR w [X]g\l — w)dw 

Jo 

= j VaR w [X]<j>(w)dw 
Jo 

= v([0,l])EX+ I TVaR w [X]{l-w)du{w) 


where 


= i/([0,l])EX + / TVaR w [X]dp(w), 
Jo 


dp(w) = (1 — w)dv(w). 


It can be shown that p is a probability measure. In fact, 


dp(w) = / u([0,w})dw 


4>(w)dw = / g'(w)dw 


1 . 


We now give some examples of interesting distortion functions and risk measures. 

Example 3.2 If w\, w 2 , w 3 , u> 4 > 0, Wi = 1, then 

g a p(x) = w^pix) + w 2 v a (x) + w 3 ^p(x) + w 4 i/; a (x), 

is a distortion function, where up, u a , 'ipp, il> a are the distortion functions of TVaR and VaR 
at confidence levels /3 and a, respectively. Then the corresponding risk measure 

Pg a p[X\ = w±TVaRp[X] + w 2 TV aR a [X] + w 3 VaRp[X] + W 4 VaR a [X\, 

is called the Glue VaR risk measure, which were initially defined by Bclles-Sampera et al. 
(2014a) (in the case = 0) and the closed-form expressions of GlueVaR for Normal, 
Log-normal, Student’s t and Generalized Pareto distributions are provided. Two new 


13 


proportional capital allocation principles based on Glue VaR risk measures are studied in 
Bclles-Sampera et al. (2014b). 

Although Glue VaR has superior mathematical properties than VaR and TVaR, how¬ 
ever, the Glue VaR risk measure may also fails to recognize the differences between two 
risks. For example, consider two risks X and Y in Example 3.1, we have computed that 
VaRo .95 [X]=VaRo .96 [X] = 100 , VaRo .95 [Y] = VaRo .96 [Y] = 100 . TVaRo .95 [X] =TVaRo .95 [Y] =300, 
TVaRo.96[X]=TVaRo.96[Y]=350. So that when a = 0.95 and (3 = 0.96, we have p 9afj [X] = 
p gafj [Y]. Thus according to p 9ap , both X and Y bear the same risk! However, the maximal 
loss for Y (1100) is more than double than for loss X (500), clearly, risk Y is more risky 
than risk X. 

Example 3.3 Let A E [0,1], define a distortion function 

g x {x) = Xg 0 (x) + (1 - A)p(x), 

where go(x) : = 1( X>Q ) and g is an arbitrary distortion function. Note that g\ can be 
rewritten as 

Jo, if x = 0, 

g\(x) = < 

| A + (1 — A )g(x), if 0 < x < 1. 

In particular, if g{x) = x, then we get the esssup-expectation convex combination distor¬ 
tion function with weight A on the essential supremum, which was introduced in Bannor 
and Scherer (2014). The corresponding risk measure 

p gx [X] = A esssup(Y) + (1 - A )E(X), 

which is a convex combination of the essential supremum of X and the ordinary expecta¬ 
tion of X w.r.t. P. 


If 


g(x) = 


[a(l -/3) + ^]x, if 0 < x < 1 — p, 

a/3 + 1 — a + a( 1 — /3)x, if 1 — p < x < 1, 

where 0<a,/3<l,0<p<l are constants, then we get 

p gx [X] = A esssup(X) + (1 - A )p g [X], 


where 

p g [X] = a( 1 - (3)E(X) + apVaR p [X] + (1 - a)TVaR p [X\. 
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As illustration, we consider the risks X and Y in Example 3.1, if p = 0.95, then 
p g [X] = p g [Y) = 50 a(3 — 250a + 300. It follows that 

p gx [X} = 500A + (1 - A)(50a/3 - 250a + 300) 


and 

p gx [Y] = 1100A + (1 - A)(50a/3 - 250a + 300). 

Taking A = a = 1 ,/3 — 0, then p gx [X] = 275 and p gx [Y] = 575. Taking A = a = /3 = 
then p gx [X] = 437.5 and p gx |Y] = 737.5. Thus the measure p gx can measure the differences 
between two risks X and Y. 


3.3 A copula-based approach 

If F is a distribution function on [0,1], then F can be used as a distortion function. 
The well-known examples are the PH transform and the dual power transform and, more 
generally, the beta transform; see Wrich and Hardy (1999) for details. Similarly, we use 
this technique to a distribution function on [0, l] 2 . We first introduce the notion of copula 
in the two-dimensional case. 

Definition 3.1. A two-dimensional copula C(u,v ) is a bivariate distribution on the 
square [0, l] 2 having uniform margins. That is a function C(u,v) : [0, l] 2 — )■ [0,1] is right- 
continuous in each variable such that (7(0, v ) = C(u , 0) = 0, C(u, 1) = u, (7(1, v) = v and 
for u 2 > ui,v 2 > t’i, 

C(u 2 , v 2 ) - C(u 2 , v\) - C(ui, v 2 ) + C(ui, vi) > 0. 

For an introduction to copula theory and some of its applications, we refer to Joe (1997), 
Denuit et al. (2006) and Nelsen (2006). 

The well-known examples of copulas are C + (u,v) = min(w, v), (7 ± (u, v) = uv and 
C~(u,v) = max(« + r — 1,0) describing, respectively, comonotone dependence, indepen¬ 
dence and counter mo notone dependence between two random variables A" and Y. The 
copula version of the Frechet-Hoeffding bounds inequality tells us 

C~(u,v ) < C(u,v ) < C + (u,v). 

Any copula has the following decomposition (cf. Yang et al (2006)) 

C(u, v ) = aC + (u , v ) + /3C ± (u : v ) + 7(7 _ (m, v ) + lG(u , u), 


15 


where a, /?, 7 , l > 0, a + [5 + 7 +1 = 1. Here G is a copula which called the indecomposable 
part. 

For a given two-dimensional copula G(-, •), define one-parameter family {g p }p e (o,i] by 
g p (u) = ( ( “' p) or c ^' u) . Clearly, for each p, is a right continuous distortion function. 
For example, 


• Sp { u ) — C = u is continuous and both convex and concave, the associated risk 
measure is EX: 

• gi-p(u) = — = min|j^,l| is continuous and concave, the corresponding 

risk measure is TVaR p ; 


g p (u) = -—= max | u+ ^ 1 ,0 J- is continuous and convex, the corresponding risk 


measure is p 9p [X] = - J Q P VaR q [X]dq. 


Conversely, if {(?p}p G (o,i] is a family of distortion functions, then, however, C(u,p ) = 
pgp(u) is not a copula in general; A sufficient condition can be found in Cherubini and 
Mulinacci (2014). 

A lot of copulas and methods to construct them can be found in the literature, for 
example, Joe (1997), Denuit et al. (2006) and Nelsen (2006). We give below the most 
common bivariate copulas and the corresponding distortion functions. 


• The Archimedean copulas: 

C 9 (u,v) = T[- 1 l(T(u) + T(u)) 

for some generator T : (0,1] —y M + with v k(l) = 0 such that T is convex. The pseudo¬ 
inverse of T is the function TH with Dorn'lb” 1 ! = [ 0 , 00 ] and Ran\l/H = [ 0 , 1 ] given 
by 


t ) = 


if 0 <t<T( 0 ), 


0, 


if T( 0 ) <t< 00 . 


If T is twice differentiable and 4/(0) = 00 , then C $ is componentwise concave if, and only 
if is concave, where 4/' is the derivative of T (see Dolati and Nezhad (2014)). Aa a 
consequence, we have 


16 








Theorem 3.1. For each v > 0, the distortion function 


g v (u) = -4^(4^) + T(u)) 
v 

is concave if, and only if -A is concave. 

We list some examples of the Archimedean copulas and the corresponding distortion 
functions: 

(a) The Clayton copula with parameter a > 0 is generated by T(f) = -( t~ a — 1) and 
takes the form 

C a (u, v ) = (u~ a + v~ a - 1)" 1/Q , a > 0. 

The limit of C a (u,v ) for a 0 and a f oo leads to independence and comonotonicity 
respectively (Nelsen, 2006). The corresponding distortion functions: 

g a , v {u) = -( u~ a + v~ a — a>0,i)6 (0,1]. 

In particular, if a = 1, we get the proportional odds distortion which is found by Cherubini 
and Mulinacci (2014): 

9i,v{U) — /h \ • 

u + u(l — u) 

Since (■^y) // = — a(l + a)t “ _1 < 0, g a , v (u ) is concave. 

(b) In case T(f) = — In ^ e e -all j we get the Frank copulas: 

- l)(e~ av - 1 ) 


C a (u,v) = — — In ( 1 + [e 
a 


The corresponding distortion functions: 


e~ a - 1 


, 0. 


1 / (e 

9a,v( U ) = - ln 1 + 

av ' 


(e~ au - l)(e~ av - 1 ) 


e~ a - 1 


Since (■^y) ,/ = —ae at , g a ,v(u) is convex if a < 0 and concave if a > 0. 

(c) In case T(t) = t -1 /" — 1 we get the Pareto survival copulas: 

C a (u, v ) = (u~ l / a + v~ 1/a - 1 ), a > 0 . 

The corresponding distortion functions: 


g a ,v(u) = - (u 1/01 + v 1/a - 


i)' 
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Since (^y)" = -(1 + < 0, g a , v (u ) is concave. 

(d) In case T(t) = In 1+ ° : ^'~ 1 - > , a G [—1,1) we get the Ali-Mikhail-Haq copulas: 

Cg(u, v) = 


uv 


1 — a(l — w)(l — v) 


The corresponding distortion functions: 

9e,v( u ) = 


u 


1 — a(l — w)(l — n) 

Since = ^rj, ge,v( u ) is convex if a G [—1, 0] and concave if a G [0,1). 

(e) In case T(f) = (— lnf) Q we get the Gumbel-Hougaard copulas: 

C a (u,v) = exp j— ((— In u)° + (— lnn) a ) 1//a |, a > 1 . 

The corresponding distortion functions: 

9a,v 0 ) = - exp (- ((— In u) a + (-lnn)") 1/Q 


The value a = 1 gives independence and the limit for a ^ oo leads to comonotonicity. 
Since 


= ---(— lnt)“ -2 ( —1 + — —— ' ^ -n ;r„2-a 

Vv h'(t) ; at \ In t 

< 0 , if 0 < t < e 2 ~ a , a > 2 , 


< 0 , if 0 < t < 1,1 < a < 2 , 
> 0 , if e 2 ~ a < t < 1 , a > 2 , 


9a,v(u) is concave if 1 < a < 2 and, if a > 2 , g a , v ( u ) is convex on (e 2 Q , 1] and concave 
on ( 0 , e 2-Q ). 

Among other copulas, which do not belong to Archimedean family, it is worth to 
mention the following three copulas, given in the bivariate case as: 


• The Farlie-Gumbel-Morgenstern copulas: 

C a (u, v) = uv + uva( 1 — u)(l — v), a G [—1,1], 
The corresponding distortion functions: 

g a , v (u) = u + ua( 1 — u)(l — v), a e [—1,1], v G [0,1], 
which is convex if a G [—1, 0] and concave if a G [0,1]. 
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The Marshall-Olkin copulas: 


C a ^(u,v) = min-fu 1 a v,uv 1 /3 }, a,/3 G [0,1]. 

Note that this copula is not symmetric for a 7 - (3. The corresponding distortion functions: 

9a,t3,v( u ) = min{u 1 _a ,wr^}, a,P € [0,1], v € (0,1], 
which is concave. In particular, g a ,o,v( u ) — 9o,p,v( u ) — u , 9i,i,v( u ) = minjl, ^}. 


The normal copulas: 


C P {u,v) = $p(<I> 

where <f> p is a bivariate normal distribution with standard normal marginal distributions 
and the correlation coefficient — 1 < p < 1 , $ _1 is the inverse function of the standard 
normal distribution. The corresponding distortion functions: 

9p,v( u ) = ^$ P ($ -1 (u),$ -1 (u))- 

4 Tail-asymptotics for VaR 


Subadditivity is an appealing property when aggregating risks in order to preserve the 
benefits of diversification. Subadditivity of two risks is not only dependent on their 
dependence structure but also on the marginal distributions. Value at risk is one of 
the most popular risk measures, but this risk measure is not always subadditive, nor 
convex, exception of clliptically distributed risks. This family consists of many symmetric 
distributions such as the multivariate normal family, the multivariate Student-f family, 
the multivariate logistic family and the multivariate exponential power family, and so on. 
A recent development in the VaR literature concerns the subadditivity in the tails (see 
Danfelsson et al (2013)) which demonstrate that VaR is subadditive in the tails of all 
fat tailed distributions, provided the tails are not super fat. However, in most practical 
models of interest the support of loss is bounded so that the maximum loss is simply 
finite. We will also show that for this class losses VaR is subadditive in the tail. We can 
illustrate the ideas here with three simple examples. In Examples 4.1 and 4.3, X and Y 
are independent, while in Example 4.2, X and Y are dependent. 
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Example 4.1 Let X and Y be i.i.d. random variables which are Bernoulli (0.02) 
distributed, i.e. P{X = 1) = 1 — P(X = 0) = 0.02. Then 

P(X + Y = 0) = P(X = 0 )P(Y = 0) = 0.98 2 = 0.9604, 


P(X + Y = 1) = P(X = 1 )P(Y = 0) + P(X = 0 )P(Y = 1) = 0.0392, 

P(X + Y = 2) = P(X = 1 )P{Y = 1) = 0.0004. 

Dhaene et al. (2006) verified that VaR is not subadditive since VaRo. 975 [V] = VaRo. 975 [V]=0 
and VaR 0 . 975 [V + Y]= 1. However, for p > 0.98, VaR p [X]=VaR p [V]=l and 


VaR p [X + V] 


1, if 0.98 < p < 0.9996, 

2, if p> 0.9996. 


Thus for p > 0.98, 


VaR p [X + Y] < VaR p [X} + VaR p [Y}. 


Example 4.2 Suppose we have losses X and Y , both dependent on the same under¬ 
lying Uniform(0,l) random variable U as follows. 


1000, if U < 0.04 
0, if U > 0.04 


Note that 


Y 


0, if U < 0.96 

1000, if U > 0.96 


P(X + Y = 0) = P(X = 0, Y = 0) = P(U > 0.04, U < 0.96) = 0.92, 
P(X + Y = 1000) = P(X = 0, Y = 1000) + P(X = 1000, Y = 0) = 0.08. 


Hardy (2006) found that VaRo.95[V]=VaR 0 .95[V]=0, VaR 0 . 95 [V + V]=1000. Thus 


VaR om [X + y] > VaR 0 . 95 [X] + VaR 0 . 95 [Y], 


However, for any a > 0.96, VaR Q [V]=VaR Q [V]=1000, VaR Q [X + Y] = 1000. Thus, 


VaR a [X + y] < Va,R a [X] + VaR a [Y). 
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Example 4.3 Let A" and Y be i.i.d. random variables which are Uniform(0,l) dis¬ 
tributed. Then 


0 , 


Fx+y(z) = < 


I~2 
2 ’ 


if z < 0, 
if 0 < z < 1, 
1 — |(2 — z) 2 , if 1 < 2 < 2, 
1 , if z> 2, 


and for p E (0,1], VaR p [X] = VaR p [Y] = p, 

V2P, 


V aR p [X + Y] = 


if pE (0, §], 


2 — a/2(i — p)i ifp e [|, 1]. 


Thus for p E 1], 


VaR a [X + Y] < VaR a [X] + VaR a [Y}. 


Generally, we have the following result. 

Theorem 4.1. If the risks Xi,X 2 , • • • , X & have finite upper endpoints, then 


lim sup 


VaRp\El iXi. 

i VaRplXi] 


< 1 . 


Proof The proof is very simple. Denote by esssup(Aj) = sup{x : P{X t < x) < 1}. 
Then esssup (Aj) < oo and P(A* < esssup (Xj)) — l,i — 1,2, ■ ■ ■ ,k, which lead to 


P 


< 'y ' esssup (X,) 


k 1=1 


i =1 


l. 


Hence 


and the result follows. 


esssup 



k 

< esssup (A,;), 

i =1 


Remark 4.1. Many distributions, such as Binomial, Uniform, have finite upper end¬ 
points; Any truncated distribution: whether it is right truncated or doubly truncated all 
have finite upper endpoints. 


Next theorem consider the random variables A 1 ,A 2 ,-- - , X& that are not necessar¬ 
ily has finite upper endpoint, we first recall the notion of (extended) regularly varying 
function: 
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Definition 4.1. A function f is called regularly varying at some point x (or x + , respec¬ 
tively) with index aeK if for all t > 0, 


lim ijA = f 

Ax f{s ) 


(or lim,j x yAi = t a , respectively). We write f E 7 Z* (f E FZ X ^, respectively). For a = 0 
we say f is slowly varying; for a = —oo rapidly varying. 


Definition 4.2. Assume that F is the distribution function of a nonnegative random. We 
say F belongs to the extended regular variation class, if there are some 0 < a < (5 < oo 


such that _ _ 

s'* < lim inf < lim sup 3^- < for all s > 1 , 

x^-oo F(x) xm-oo F(x) 


or equivalently 


s " < lim inf -=■—- < lim sup — ' ^ < s 13 , for all 0 < s < 1 . 

x^oo F(x) moo F(x) 


We write F E ERV(—a, —f3). 


A standard reference to the topic of (extended) regular variation is Bingham et al. 
(1987) while main results are summarized by Embrechts et ah (1997). 


Theorem 4.2. We assume that Xi,X 2 ,--- , X/,. have the same absolutely continuous 
marginal distributions F with infinite upper endpoint. 

(1) If 



Um P (Ai Xi>z ) K 1 
- P(X L > f) < ’ 

then 

'” 1 Eli VaR, [X,] 

(2) If 

iim P( EL,V>.-) =1 

P(X 1 > f) 

then 

.■VaR, EhiXJ , 

P - 1 E?=i VaRplXi] 


(4.1) 


(4.2) 
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(3) If 

lim 


2—>-oo 

then 

lim 


P (Hi=l X i > z ) 

nx i > i) 

VaRplEl iXj_ 
^ Ef=i VaRpiXi] 


> 1 , 


> 1 . 


Proof We prove (1) only since the other cases follow immediately in the same way. 
Because all the marginal distributions are absolutely continuous, so we have for any 
pe ( 0 , 1 ), 


P(X i > VaRpiXi]) =P[j2 X i> VaR v 


j= i 


Y. X < 


This, together with (4.1), implies that 


P(Xi > VaRp[Xi]) 


lim 


1=1 


< 1 . 


= 1 — p. 


(4.3) 


The absolute continuity of F implies that F is continuous and strictly monotone decreas¬ 
ing. Then from (4.3) we have 

lim Va ^ ] >1, 

paflpEh, A'] 

which is (4.2). This completes the proof. 


Example 4.4 Suppose that each X^ is regularly varying with index —a < 0. When 
the X t are mutually independent, it follows from ( Feller 1971, p. 279) that 


Thus we get 


lim 

s —^OO 


P(TLXi>8) 

P(x i > f) 


k_ 

k a ' 


lim VaR ,Eti A.] 

Eli VaR,[X,] 


< 1, 

if a > 1, 

= 1, 

if a — 1, 

> 1, 

if a < 1. 


Suppose that the X t are commonotonic, i.e. P(Xi = • • • = X/.) = 1, then 


lim 

s—> oo 


P(Ill 1 Xi>8) 
P(X 1 > f) 


1 . 
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So that in the case a = 1 the result for the independent and the commonotonic case are 
the same. 

The following result generalizes Theorem 10 in Jang and Jho (2007) in which all Y^s 
are assumed identically distributed. 

Theorem 4.3. Suppose Yi, • • • , are nonnegative random variables (but not necessarily 
independent or identically distributed.) If Y\ has distribution F satisfying 1 — F{x) = 
x~ a L(x),a > 0,x > 0 ? where L 6 77 “ is slowly varying at infinity. If —» c, and 


P(Yi>x,Yj>x) 

F{x) 


— y 0 , i j, as x oo, i, j = 1 , 2 , • • • , k, then 



Proof It follows from Lemma 2.1 in Davis and Resnick (1996) that 



This leads to 



(4.4) 


Because 


P(Yi > VaR p [Y i]) = P K Y t > VaR p ^ Y t 



i= 1 


k 


l — p. 


Thus from (4.4) that 


lim P(r, > = 

P(Y X > VaR ,,Ell yj) tr ■ 


which is equivalent to 


lim 


P{Y i > VaR p [Yi]) 


1 . 



This implies that 


VaR p [Yi] 


1 , 
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rsj 


since F is continuous and strictly monotone decreasing. Note that c\ = 1, cJ a VaR p \Yi\ 
VaR p [Yi\ (as p —> 1) and 



completing the proof. 

Remark 4.2. The above result is obtained by Embrechts et al. (2009) for identically 
distributed and Archimedean copula dependent ’s. However, our result can not obtained 
from their’s due to the following fact: The famous Farlie-Gumbel-Morgenstern family, 
does not belong to Archimedean family, which has the form 


F(x,y) = F 1 (x)F 2 (y)(l + aF 1 (x)F 2 (y)) 


where F X ,F 2 are two distributions and a G [—1,1] is a constant. When F 1 = F 2 , it 


satisfying —» 0 as x —> oo. 


In the next theorem we consider the extended regularly varying instead of regularly 
varying. 

Theorem 4.4. Suppose Yi, ■ ■ ■ , Y k are nonnegative random variables with the common 


identical distribution function F. If F e ERV(—a,—/3 ) and 
x —$■ oo, i,j = l,2,---,k, then 

(1) If (3 < 1 , 


P(Yj>x,Yj>x) 

F(x) 


-t 0 ,i yi j, as 



(1) If a > 1, 



(1) If a = 13 = 1, 


E*_i VaR p \Y,] 



Proof It follows from Lemma 2.2 in Zhang et al. (2009) that 



—>■ k, as x —^ oo. 
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This leads to 


lim p (Eti Xl > yafl pEti UP 
P(Y i > VaBjlEl., VJ) 


= k, 


f-i p(y, > VoflpEli *]) 

,. FfrjyoiipEt!^]) 

lim sup--- - — 1 


= k. 


from which and using the same argument as that in the proof of Theorem 4.3 leads to 

]im P{Y, > Va iyy]) 

If f3 < 1, then 

~F(]c~h/nP \S^ k VA\ 


This and (4.5) imply that 


^ (k~?)-PF(VaR p Eli^D 


hm sup = -- 


F(VaR p [Yi\) 


It follows that 


Thus 


LJ.X kJ L *T _ 1 _ ; 

p - 1 Fik—pVaR^YH^Y,]) 

,. k^VaRpiY,] ^ i 

hmsup- t-< 1. 

p-i VaR.lEt^ ~ 

y kV a R, P [Y]\ 

hm sup- 


> 1. 


VaRplYlt^Yt] 


< 1. 


Similarly, if a > 1, 


lim inf 


and hence 


k^VaRplYi] 

EC VaiSpEl,! Y,] 

limmf WaR f l] 

VaRJSXM 

If a — /3 = 1, then by (4.6) and (4.7) one has 

, kVaRJYi} 

hm-- 

This ends the proof of Theorem 4.4. 


> 1 


> 1. 


= 1. 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


To give applications of our results we employ extreme value theory techniques. A 
distribution function F (or the rv X) is said to belong to the Maximum Domain of 
Attraction (MDA) of the extreme value distribution FI if there exist constants c n > 
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0, d n G M such that c“ 1 ({max(Xi, • • • , X n j — d n ) - 4 - H. We write X G MDA(H) or F G 
MDA(H). According to the Fisher-Tippett theorem (see Theorem 3.2.3 in Embrechts et 
al. (1997)) H belongs to one of the three standard extreme value distributions: 

Frechet type : $ a (x) 

Weibull type : T a (x) 

Gurnbel type : A(x) 

Let xp denote the right-endpoint of the support of F: Xp — inf (x : F(x) = 1}. Then we 
have the following results (see Embrechts et al. (1997), PP. 132-157). 

• Frechet case: For some a > 0, F G MDA{fb a ) 77°° a . 

Examples are Pareto, Cauchy, Burr, Loggamma and Stable with index j3 < 2. 

• Weibull case: For some a > 0, F G MDA(f^ a ) o xf < oo, F(xf — 1/x) G 77°° Q . 
Examples are Uniform and Beta distribution. 


0, if x < 0, 

a > 0. 

exp{— x “}, if x > 0, 

exp{ — (— x) a }, if x < 0, 

1, if x > 0, 
exp{— e~ x \, i6l 


a > 0 . 


• Gurnbel case: F G MDA( A a ) xf < cxo and there exists a positive measurable 

function a such that for t G R 


lim 

x'lxp 


F(x + ta(x )) 
F(x) 


= e 


—t 


(4.9) 


Examples are Exponential-like, Weibull-like, Gamma, Normal, Lognormal, Benktander- 
type-I and Benktander-type-II. 


Remark 4.3. (1). For a > 0, if Xi, X 2 , ■ ■ ■ ,Xk G MDA(f$> a ), in view of Weibull case 
above they are all have finite supports. It follows from Theorem f.l, VaR p is subadditive 
forp is sufficiently close to 1. 


(2). For a > 0, if X\ , X 2 , ■ ■ ■ , X k G MDA(& a ) and are identically distributed, 
(-X 1 ,—X 2 ,--- ,—X k ) has an Archimedean copula with generator if, which is regularly 
varying at 0 with index —/3 < 0. We apply (2.2) in Alink et al. (2004) an d Definition 4-1 
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to obtain 


lim 

z—> oo 


n T,l 1 x i >z) 

n* i > f) 


P(Eti E > *) P(Xi > 0 

p(x i > f) 

P(E > z) 


lim 

2^-00 P(A 7 ’l > ZJ 


qk(/3,cn) lim 


lim 

Z^CO 


P(Xi > f) 
P{X 1 > z{q k {P,a))~ 1 l a ) 
P(X i > f) 


{ <1, if a > 1, 

= 1, if a = 1, 

>1, if a < 1, 


where in the last step we have used Lemma 3.1(d) in Embrechts et al. (2009) which states 
that 

min{fc a , k} < g fc (/3, a) < ma x{k a , k}. 


This, together with Theorem 4.2 we recover the result Theorem 2.5 in Embrechts et 
al. (2009). 


(3). If Xi, X 2 , • • • , Xk € MDA(A a ) have common distribution F, (—Xi, —X 2l ■ ■ ■ , — Xk) 
has an Archimedean copula with generator -0, which is regularly varying at 0 with index 
—/3 < 0. We apply (2.6) in Alink et al. (2004) to obtain 


lim 

z—too 


p (Et, v > z) 

P(X i > |) 



where 



The constant qf(/3) < e* is strictly increasing in /? with 


lim Qk (/3) = 0, lim q°(p) = e K 
/3-v 0 /3->oo 

For more details, see Alink et al. (2004) for the case k = 2 and Chen et al. (2012) for 
general case. Thus by Theorem 4.2, 


|im v«RAT.U x ,] 

»-■ Ell VaR, [.V] 


< 1. 


In particular, when a —> oo, 










Remark 4.4. Note that VaR p [X] is a left-continuous nondecreasing function having 
VaR 0 [X] as the essential infimum of X and VaR,i [X] as the essential supermum of X. 
Thus under the assumptions of Theorem f.l or Theorem 4.2, if p close to 1, we have 

VaR p [Xi + X 2 ] < VaR p [Xi] + VaR p [X 2 }, 

which, together with the positive homogeneity ofXaR p [X], implies that, if p close to 1, the 
convexity is holds: 

VaRp [AX + (1 - A )Y] < AVaR p [X] + (1 - A)VaR p [V], 0 < A < 1. 


From above analysis we see that, although, in general the VaR risk measure lack of 
subadditivity and convexity. However, one should not too worries about violations of 
subadditivity for risk management applications relying on VaR, since in most practical 
circumstances it is subadditive, at least is subadditive in the tail, and the failure to be 
subadditive in a few situations is not sufficiently important to reject the VaR risk measure. 


5 Tail-subadditivity for distortion risk measures 


The tail-subadditivity property for Glue VaR risk measures were initially defined by Belles- 
Sampera et al. (2014a) and the milder condition of subadditivity in the tail region is 
investigated. Furthermore, they verified that a Glue VaR risk measure is tail-subadditive 
if its associated distortion function k h ^ 2 (u) is concave in [0,1 — a), where parameters a 
is confidence level and /3 is an extra confidence level such that 0 < a < j3 < 1 and, 

{ j^u, if 0 < u < 1 — /?, 

h\ + h p_^f (u — 1 + (3), ifl — (3<u<l — a, 

1, if 1 — a < u < 1, 

where hi and h -2 are two distorted survival probabilities at levels 1 — /3 and 1 — a, respec¬ 
tively. Here 0 < h\ < /12 < 1. We note, however, from their proof to Theorem 6.1 that 
the result will hold for any distortion function that is concave in [0,1 — a), not restricted 
to k h ^ 2 {u). In this section we state the result and give an alternative proof. As in 
Belles-Sampera et al. (2014a), for a given confidence level a, the tail region of a random 
variable Z is defined as Q a ,z = {w\Z{w) > s Q } ^ where s a = i\\t{z\Fz{z) < 1 — a} is 
the ct-quantile. For simplicity, we use the notation Sz(z ) := Fz(z). 
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Theorem 5.1. For a confidence level a € [0,1] and two risks X,Y defined on the same 
probability space. If Q a ,x fl Q a ,Y H Q a ,x+Y 0 and g is a concave distortion function in 
[0,1 — a), then the distortion risk measure p g is tail-subadditive. That is 

pO poo 

/ [s(5x + y(z)) - 1 ]dz + / g(S x+ Y(z))dz 


' 0Am o 


< 


' 0 Vm Q 
0 

'OA m 0 


[g{S x (z )) - 1 ]dz + / g(Sx{z))dz 


' 0 Vm„ 


+ / \g(S Y (z))-l]dz+ g(S Y (z))dz , 


' OA m 0 


' 0Vm Q 


where m a = sup{s a (X), s a (F), s Q (X + F)}. 

Proof Without loss of the generality, we assume that the risks X and Y are nonneg¬ 
ative, so that m a = 5^! y (1 — a) > 0. It follows that 


g{S x +Y{z))dz = / g{S x +Y{z))dz 

eOO 


dx 


’S x+Y {l-a) J[ 0 ,S x +y(x)} 


dg(g) 


^ X+Y (■*■ S') 


- x+yv 

= / / eh 

^[0,Sx+y(S x l,. y (l-a))) ^Sx+y( 1_Q ) 

J [0,Sx+y(S' J¥ ^y(l-Q;))) 

where in the third step we have used the Fubini’s theorem to change the order of integra¬ 
tion. As above, we have 

r r i 

r>-l 


'[0,Sx+y(S x + y (l-a))) 


F x+Y (l-q)dg(q) = 


Fx+v^d'O- - q) d q 


ll-S x +Y(S x+Y (l-a)) 

/ TVaR x +Y(w)dnx+Y(w ), 

' o 


where 


d/i X +y(w) = (1 - w)di/x+y(^), ^x+y([0,g]) = l( 1 -s jr+y( s^ y (i-a)) I i](9)^( 1 ~ <?)• 

Finally, we get 


g(S x +Y{z))dz= j TVaR x+Y {w)dg x +Y{w) - S x l +Y (l-a)g(l-a). (5.1) 

Jo 


Similarly, 


POO pi 

g(S x (z))dz= TVaR x {w)dnx{w) - S X \ Y {1-a)g(l-a), (5.2) 

7l a J 0 
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and 


where 


and 


g(Sy(z))dz — / TV aRY(w)diiY(w) — S x+Y (1 — a)g(l — a), (5.3) 


df*x(w) = (1 - w)du x (w), vx([0, q\) = l (1 _ Sjf(s -^ (1 _ a)))1] (<?)£/'(! - q), 


dfj, Y (w) = (1 -w)du Y (w), z/y([0 ,q]) = l {1 _s Y{s -i +Y(1 - a))A] (q)g\l-q). 

By the subadditivity of TVaR and note that vx+y ([0, q]) < vx([ 0, <?]), ^([0, q]), we obtain 

[ TVaR x+Y {w)dnx+Y(w) < [ TVaR w (X)dnx+Y(w) 

Jo Jo 

+ / TVaR w (X)dnx+Y(w) 

Jo 

< / TVaR x {w)dnx{w) 


+ / TVaRY(w)d/iY(w), 
Jo 

this, together with (5.1)-(5.3), implies that 

poo poo poo 

/ g(S x +Y{z))dz < / g(Sx(z))dz + / g(S Y (z))dz, 
J m a J m a J m a 

as desired. 


Remark 5.1. Consider the distortion functions associated with the Gumhel-Hougaard 
copulas (cf. Section 3.3): 

g a ,v{u) = ^ exp j- ((— lntt) a + (— lnu)“) 1/Q j . 

If a > 2, then g a ,v(u) is concave on (0,e 2 ~") and convex on (e 2-Q , 1]. Thus the distortion 
risk measure p ga v ( u ) is tail-subadditive. 
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